Conformal Transformation and the Complex Potentials in 

Ideal Hydrodynamics 


M.Cattani 

Institute of Physics of the University of Sao Paulo - Brasil 

(mcattani@if.usp.br ) 


Abstract. 

This article was written to graduate and postgraduate students of 
physics and engineering. We show how complex functions, conformal 
transformation and complex potentials are used in hydrodynamics of ideal 
fluids. We have also briefly analyzed complex potentials in electrostatic. 
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(I)Introduction. 

Complex functions, conformal transformation and complex 
potentials are very useful mathematical tools that can be used to describe 
many physical phenomenon like, Hydrodynamic and Heat flows, Elasticity, 
Gravitational Field, Electrostatics, Magnetism and Electric Current flows. 
Our paper is a didactical one written to students of physics and engineering 
showing how hydrodynamics problems can be studied using the theory of 
functions of complex variable. There are excellent text books on this 
subject but, to be concise, we will consider only the Irving and Mullineux 
book. [1] In Section 1 are defined Analytic Functions. In Section 2 we 
present Conformal Transformations that connect complex planes. In 
Section 3 are defined analytic functions, named Complex Potentials. In 
Section 4 is seen how some Hydrodynamic Flows of ideal fluids can be 
studied using the analytic formalism defined above. In Section 5 is shown 
the Uniform Flow around a Plate. In Section 6 is seen the Flow Around an 
Airfoil. Finally, in Section 7 we briefly show how 2-dimensional 
Electrostatic Problems can be studied with the approach seen above. 

(1) Function of a complex variable. 

Let us consider a two 2-dimensional complex plane with variable 
z = x + iy and an analytic function co(z), expect at a finite number of 
singularities^ 1] 

co(z)= f(z) — u(x,y) + i v(x,y) (1.1), 
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( 1 . 2 ), 


where u and v obey the Cauchy-Riemann equations, [1] that is, 

du/dx = dw/d y and du/dy = - owl ox 

and the 2-dim Laplace equations, 

o 2 u/ox 2 + o 2 u/oy 2 = 0 and d 2 w/dx 2 + d 2 v/d y 2 =0 (1.3), 
showing that u(x,y) and v(x,y) are harmonic functions, and co(z): 

dco(z)/dz = (du/dx) + i (owlox) - (cw/dy) - i (du/dy) (1.4). 

(2) Conformal transformation^ 1] 

The equation (1.1), w = F(z) connects two 2-dimensional complex 
planes with variables z = x + iy and co(z) = F(z) = u + i v, respectively. It 
can be interpreted as mapping given locus in the z-plane into the locus in 
the co-plane. It is possible to find a wide class of such transformation in 
which complicated loci in the z-plane are mapped into fairly simple ones in 
the e-plane. It is by this method many problems of physics are solved. 
Differentiation of the function s = F(z) gives 

dco/dz = F'(z) = x exp(ia) (2.1), 

where % = IF'(z)l is supposed neither zero nor infinite. Equating the moduli 
and phases of each side of (2) we have 

ldco| = x |dz| and ph( dco) = a + ph( dz) (2.2), 

where ph(...) means "phase" or "argument" of (...). Let us suppose that in 
the z-plane we have a curve (locus) T and that in the co-plane the 
corresponding figure (locus) y. Equations (2.2) mean that: 

(A) An infinitesimal arc of y is % times the corresponding arc of Y, that is, 
the y arc is magnified by % factor y and that the arc ds is turned through an 
angle a in the positive sense; 

(B) The phase angles under the transformation are preserved in magnitude 
and sense. 

To see this we note that in the z-plane near Pj we have dz 2 /dz! = 
F'(Zi), that is, 0j= ph( dz 2 ) - ph(dzd) = ph[F'(zd}] Pl where 9i is the 
difference of angles of the respective dz's. Similarly, in co-plane dco 2 /dcoi = 
E'(zi), that is, 0 2 = ph( dco 2 ) - ph( coO = ph[ F'(zi)] P i. where 0 2 is the 
difference of angles of the respective dco's. 

This complex transformation, at all nonsingular points, is said 
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to be conformal: under it infinitesimal figures in z-plane transform into 
similar, but re-oriented, figures in the co-plane. 

Note that the coordinate transformation between two manifolds is 
valid only when it is possible to obtain the inverse transformation. That is, 
when the Jacobian of the transformation does not vanish or it is not infinite. 
That is, if at a point z = zi ("critical point") of the curve T we have 

IF'(zi)l = 0 or oo (2.3). 

the transformation ceases to be conformal and no single-valued inverse of 
co = F(z) exist. As co = f(z) is assumed to be regular, u and v satisfy the 
Cauchy-Riemann equations. In this way, we verify that the Jacobian of the 
transformation is given by, 

J(x,y,u,v) = (5u/5x) 2 + (dv/dx) 2 = IF'(z)l 2 (2.4). 

When IF'(z)l f 0 or oo, the transformation is conformal. 

As a consequence of the item (B), orthogonal families of curves 
u(x,y) = constant and v(x,y) = constant in the z-plane are also orthogonal 
families of curves in the co-plane. 

Note that if ®(x,y) = r|(x,y) + i £,(x,y) is a solution of the Laplace 
equation, and if under the conformal transformation F(z) it becomes 0(rj,^) 
then ®(r|,£,) also satisfy the Laplace equation, that is, 

0 2 O(r|^)/dr| 2 + 0O 2 (r|^)/^ 2 = 0 . 

Example (2.1) 

Straight lines x = constant and y = constant in the z plane and F(z) = cos(z). 

So, writing 

cos(z) = (e lz - e' lz )/2 = (e lx ' y - e' lx + y )/2 = cosh(y)cos(x) + i sinh(y)sin(x), 
that is, cos(z) = u(x,y) + iv(x,y), where 

u(x,y) = cosh(y)cos(x) and v(x,y) = sinh(y)sin(x). (2.5). 

From Eqs.(2.5) we get the following curves in the co-plane: 

u 2 /cosh 2 (y) + v 2 /sinh 2 (y) =1 and u 2 /cos 2 (x) - v 2 /sin 2 (x) = 1. 
That is, curves with x = constant in the z-plane become hyperboles in the 
co-plane with focus in u = -1 and u = 1; curves in z-plane with y = constant 
become ellipses with the same focus (Figure 1). We know from analytical 
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geometry that these conics with the same focus form a system of 
orthogonal curves. 



Figure 1. Straight lines x = const and y = const in the z-plane and conformal figures in 

the co-plane when F(z) = cos(z). 


(3) Complex Potentials. 

Let us define Complex Potential Q(z) by the analytic function 

Q(z) = F(z) = cj)(x,y) + i \)/(x,y) (3.1), 

where (J)(x,y) and \|/(x,y) are real harmonic functions that obey the relations 
seen in Sections land 2. That is, dtydx. = dy/dy, d\\f/dx = - dfy/dy and 
A(x,y) = A\|/(x,y) = 0, where A = Laplacian operator. If these functions 
represent orthogonal families of curves <J)(x,y) = constant and \|/(x,y) = 
constant in the z-plane they continue to be orthogonal families of curves in 
the Q-plane defined by (3.1). Due to these peculiar properties very useful 
information are obtained interpreting physically these functions in various 
branches of applied mathematics in Physics. This is shown in the Table I. 



4>( X 'V) 


'/>(*. y) 

H yd rod y n am i cs 

Velocity potential 


Stream function 

Heat flow 

Isothermals 


Heat flow lines 

Elasticity 

Strain function 


Stress lines 

Gravitational field 

Gravitational potential 

1 


Electrostatics 

1 

Lines of force 

Magnetism 

Current flow 

f Potential 

Lines of flow 



4 






























(4) Hydrodynamics. 

Some of the most interesting applications of the this theory are to 
problems of ideal fluids hydrodynamics, and so the tendency will be to 
interpret the complex potential applicable to this subject. 

Indeed, in hydrodynamics a stream line is defined by a stream 
function \|/(x,y) that obeys the conditions [2] 

V x = dy/dy and V y = - d\\i/dx, (4.1) 

where V x and V y are the components of the velocity of the stream parallel 
to the x and y axes, respectively. Along a stream line we have 

d\|/ = (8\\f/6x)dx + (8\\i/8y) = 0. (4.2) 

As, for an incompressible fluid div(V) = 0we see that V can also be written 
in terms of a velocity-potential function c()(x,y): 

V x = dtydy and V y = - dtydx. (4.3). 

From (4.1)-(4.3) we also verify that 

A())(x,y) = A\|/(x,y) = 0 and grad((() )«grad(\|/) = 0 (4.4), 

that is, § and \|/ are harmonic functions and that the stream lines 
(represented by a continuous lines) are orthogonal to the velocity-potential 
lines (represented by a dotted lines). Thus, in hydrodynamics the potential 
function Q(z) is defined by 

Q(z) — f(z) — (|)(x,y) + i v|/(x,y) (4.5). 

In the general case when we have uniform flows very far from 
immersed body the potential Q(z ) can be written as [1,2] 

Q(z) = Az + B log(z) + a n /z n (3.2), 

where A, B, a n are complex numbers. Some simple cases of flow are 
discussed in what follows. 

(4.1) H(z) = Az^ Uniform stream 

Assuming the stream to be uniform, i.e. V x and V y to be constants 

dn/dz = dfy/dx + i dfy/dy = V x - i Y y = constant 
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showing that Q(z) is a linear function of z. That is, 

V 0 z exp(ia) = V 0 (cos a + i sin a)(x + i y), (4.7) 

where V G and a are real constants. Then, 

(p = V G (x cos a - y sin a) 

\|/ = V 0 (y cos a + x sin a) and 

V x = dfy/dx = V G cos a and V y = dfy/dy = -V 0 sin a . 

So, we have an uniform flow with velocity V 0 inclined at an angle a 
relative to the x-axis. 

(4.2) n(z) = m log(z) —> Point source (m > 0) and point sink (m < 0). 

m is a real constant and z = r exp(i0). Then, 

(j) = m log r and \|/ = m0 (4.8). 

The curves when \|/ = 0 = const are the streamlines , i.e. they constitute a 
family of straight lines passing through the origin; the equipotentials ()) = r 
= const are a family of dotted concentric circles, centre at the origin (Fig.2). 



Figure 2. Point source when m > 0 and point sink when m < 0. 

Since the stream function \|/ measures the rate of flow across any 
curve, the flow O y across any curve y encircling the origin is 

® y =J y d\|/ = J 0 2lt md0 = 27im. (4.9) 

If m > 0 fluid is emanating from the origin at a rate equal to 27im and we 
have a point source at the origin. If m < 0 we have a point sink. In either 
case m is called strength. 

The circulation around the origin by curve y is given by 

C y = J y d(J) = J y md(log r) = 0 (4.10) 
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since the value of log r returns to its original value on making a complete 
circuit of y. 

A point source or sink at a point z = a is given by 
fl(z) = m log(z -a), 
where m > 0 or m < 0, respectively. 

(4.3) H(z) = - i k log(z) —> Point vortex. 

k real and z = r exp(i0). 

So, 

4> = K0 and \|/ — - k log(r) (4.11); 

the streamlines \|/ = const are concentric circles, centre the origin and the 
equipotential lines (p are dotted lines passing through the origin (Figure 3). 



Figure 3.Point vortex with intensity k counter clock horary. 

As the stream function \|/ measures the rate of flow across any 
flow ® y across any curve y encircling the origin is 

® y = J* y d\|/ = - j y K d(log(r)) = 0 

and the circulation 

C y = I y d(|) = j 0 2lt Kd0 = 27TK. 

A point vortex at z = ± a is given by H(z) = - i k log(z ± a ). 

(4.4) Superposition of flows. 

Many examples of superposition of flows and figures can be seen, 
for instance, in reference [2]. Here only a few cases will be analyzed. 


curve, the 

(4.12) 

(4.13) 
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(4.4.1) Source + sink. 

Source and sink at z = ± a are represented by the complex potential 
n(z) = nx(z) + 0 2 (z) = nix log(z - a) + m 2 log(z + a), 

where up and m 2 can be positive or negative. In Figure 4 is shown a sink at 
x = - a with intensity m ] = - m and a sink 2 at x = a with intensity rn 2 = m. 



Figure 4. Sink at x = -a and source at x = a. 

We have a "dipole" if nt| = m = - nt 2 or a "doublet" when m, = m 2 . 

(4.4.2) Dipole. 

For a dipole we have: 

H(z) = m log[(z - a)/(z + a)] , 

Then, 

§ = m logl(z-a)/(z+a)l and \|/ = m ph[( z-a)/z+a). 

When sink and source are very far apart, Izl —> go, a —► 0 and m —> oo 
(very large source ) we verify that, putting p = 2am: 

n(z)=lim{m log[(z - a)/(z + a)]}= lim{m log[(l- a/z)/(l+a/z)]} ~ -p/z (4.14) 

In Electrostatic the point sink or source are replaced by line charges 
q = charge per unit of length, and strengths m are replaced by m = 2q/K , 
where K is the dielectric constant of the medium: [1] 

Select = (2q/K) logl(z-a)/(z+a)l and \|/ e i ec t = m ph[( z-a)/z+a). 

(4.4.3) Sink + uniform stream. 

If the flow is parallel to the x-axe that is, a = 0 (See Section 4.1): 
H(z) = -V z - Q log(z), 

so, 

(J)(x,y) = -Vizi - m loglzl = -V x - Q log[x 2 + y 2 ] 



and 


\|/(x,y) = -V ph{z) - Q phiz) = -Vy - Q tan _1 (y/x). 

The streamlines that are given by [2] 

\|/(x,y) = const = -V y - Q tan _1 (y/x), seen in Figure 5, 



Fig. 1-63 


Q (z)= — ^ln z-Vz 


Figure 5. Superposition of sink with intensity Q and uniform stream with velocity V. 

(4.4.4)Uniform flow + vortex + dipole. 

The complex potential £2(z) for the combination of an uniform flow 
with velocity V G in the negative x direction, point vortex and a dipole both 
situated at the origin x = a = 0 is given by, from Eqs.(4.7),(4.9) and (4.11), 

f2(z) = -V 0 z - i k log(z) - p/z (4.15). 

If z = r exp(i 9) the streamlines are given by \|/ = Im[Q(z)] : 

\|/(r,0) = const = -V 0 r sin 9 - k log r + p sin 9/ r (4.16). 

Let us define a number a given by p = V 0 a 2 substituting this value in 
(4.16) and putting const = -k log a, we have 

V 0 r sin 9 (r - a /r) + k log r = k log a (4.17). 

In this way, it is immediate that if r = Izl = a, then Eq.(4.17) is 
identically true for all values of 9. Detailed stream lines will not be shown 
here. So, the stream line Izl = a is a circle centre the origin radius a (Figure 
6); that is, one streamline is replaced by a rigid circular boundary, which 
can be taken to represent, for instance, the cross section of an infinite 
cylinder situated in the fluid. So, Q(z) is the complex potential for a flow 
past an infinite cylinder situated in an uniform stream when a circulation 
around a vortex, given by Eq.(4.13), 27 tk is present. The potential cp(x,y) is 
given by [1] 
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(p(x,y)= Re[Q(z)] = - V 0 x + k tan _1 (y/x) - p x/(x 2 +y 2 ) (4.18). 


With the fluid velocities, V x = dtydy and V y = - dtydx, one can estimate the 
force F x and F y (Figure 6) due to an ideal fluid on the cylinder taking into 
account the Bernouilli's equation. [L2] We have F x = 0 and F y = 27ipKV 0 
which is the lifting force called known as the Kutta-Joukowsky law. [2] 


y 



Figure 6. The stream lines are omitted. We have shown only the lifting force exerted by 
the fluid on the cylinder estimated using the Kutta-Joukowsky law) 1 ' 21 


(5)Uniform flow around a plate. 

Let us see the application of the above theory in the case of a plate 
with length d immersed in a uniform flow in the z-plane (Figure 7). In 
this figure is also shown, in the 8-plane, the flow around a circumference 
with radius a. 



/ / / / /Ja. 

Voo=V=(Vcos (X,Vsen ex) 


/ / / / / £«. 
v’oo”V={Vcos CX,V$«n (X) 


Figure 7. Illustration of the method used to estimate the flow around a plate using the 
conformal transformation. [2] 
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As this flow is extremely difficult to obtain in the z-plane using 
directly the hydrodynamic Euler's equation it will determined taking into 
account a conformal transformation, supposing that the flow conditions 
very far from of the contours in the s and z-planes are the same. This is 
done in two steps: 

(I) First, using Euler's equation we determine in the e-plane the flux Q(e) 
around a circumference with radius a = 1. This is would be given by [2] 

O(e) = C(e e' ia + e ia /e) (5.1), 

where C is an adjustable parameter. 

(II) Once known Q(s) we get the flux Q(z) in the z-plane using a conformal 
transformation 121 

£ —* z = F(s) = (d/2)(s + 1/8) (5.2), 

which transforms the circumference with unitary radius in a plate with 
length d when s varies along the circumference. 

Putting Q* = Q[e(z)] we have 

dO*/dz = (dn(s)/ds) (de/dz) = A(e' ia - e ia /s 2 )(d8/dz) (5.3). 

2 1/2 

Using (5.2) we obtain s(z) = z/d ± [(z/d)~ - 1] . Since we must consider 
only points outside the plate (|e| > 1) it is enough to take only the signal + 
of the square root. In this way, de/dz = 2e /(e -l)d. In this way, we get 

dO*(z)/dz = [2C/d]{(s 2 e' itx - e ia )/(e 2 -l)} (5.4). 

Since for £ »l(very far from the origin) we must have [2] 
dO*(z)/dz = v x (z) - i v y (z) ~ Ve' m we verify that (2C/d) =Y. Consequently, 

dO*(z)/dz = v x (z) - i v y (z) = V {(e 2 e' ia - e ia )/(e 2 -l)} (5.5). 

Using Eq.(5.5) the field velocity can be completely determined: 

v x (z) = V cos(a)/[z 2 - d 2 ] 1/2 and v y (z) = V z sin(a)/[z 2 - d 2 ] 1/2 (5.6), 

showing the "critical points" x. = - d and x+ = d, where the conformal 
transformation is invalid. This velocity flux is seen in Figure 8. 


li 




(a) (b) 

Figure 8. Illustration of the method used to estimate the stream lines (a) around a plate 
using the conformal transformation. [2] In Figure (b) is shown the force of the fluid on 
the plate that are estimated in reference [2]. 

(6)Flow around airfoil submitted to an external uniform flux. 

In the general case, [2] the flow around an airfoil submitted to an 
uniform flux can be obtained using a conformal transformation slightly 
different from that seen in Section 5. It transforms circumferences in aisles 
profiles of birds or airplanes. In the z-plane with coordinates (x,y) and 
origin O (Figure 9), a circular cylinder with radius a is immersed in an 
uniform flow with velocity V. The circle with radius a is centered at the 
point Oi taken as the origin of an auxiliary coordinate system (xi,yi). 



Figure 9. Illustration of conformal transformation to analyze flow around airfoils. 

(6.a) Flow with no vortices around the cylinder. 

In this case, the airfoil profiles in e-plane from circumferences in 
z-plane are usually obtained using conformal transformations, [2] 

e — F(zj) — aj/zj + a 2 /zf + a 3 /zi +... 

In the particular case of A = 1, B = 0, ai = a 2 and a 2 = a 3 = ...= a n = 0 
(n >1) we have the Joukowsky transform, [2] that is, e^) = Z\ + a 2 /z 1? taken 
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into account in Section 5. So, we can show [2] that in the z r plane we have 
0(z!) = -V(Zi + a 2 /Zi) around the cylinder. Thus, the flow around the airfoil 
in the e-plane can be determined [2] taking into account that e = + a 2 /z!. 

(6.b) Flow with vortices around the cylinder. 

If there is a vortices with intensity T around the cylinder [2] (see 
Section 4.3) we get Q(zi) = -V(zi + a 2 /zi) - i (I727r)ln zi. [2] 

The resulting stream lines for the cases (6.a) and (6.b) are shown, 
respectively, in Figures (9.a) and (9.b). Note that in Fig.9 the incident flux 
V and the vortex orientation T are inverted, that is, V —> - V and T —> - T. 




(9.b) 


Figures (9.a) and (9.b). Streamlines around airfoils when T = 0 (a) and T ± 0 ( b). 

We can shown [2] that the resulting forces ("drag forces") on the 
airfoils are equal to F a = 0 and F b = F L = p V T ("lift force"), where p is the 
fluid density. 


(7) Two-dimensional electrostatic problems. 

The method of complex variable potential combined with conformal 
transformation is also an especially powerful method of solving two- 
dimensional electrostatic problems. [3] In this case the complex potential 
Q(z), with z = x + iy, is written, following Panofsky& Phillips [3] as 

Q(z) = ®(x,y) + i y(x,y) (7.1), 

where ®(x,y) is the electrostatic potential and \|/(x,y) are the "stream 
functions" connected with electric field E = - grad(O). The analytic 
function f2(z) satisfy the Cauchy-Riemann relations 

ch| f/dx = - cXD/dy and 3\\i/8y = o®/ox (7.2), 


and, consequently, grad(®)* grad(\|/) = 0 and A®(x,y) = A\|/(x,y) = 0 where 
A = Laplacian operator. 
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